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The features of hot and dense gas of quarks which are considered as the quasi-particles of the 
model Hamiltonian with four-fermion interaction are studied. Being adapted to the Nambu-Jona- 
Lasinio model this approach allows us to accommodate a phase transition similar to the nuclear 
liquid-gas one at the proper scale and to argue an existence of the mixed (inhomogeneous) phase of 
vacuum and normal baryonic matter as a plausible scenario of chiral symmetry (partial) restoration. 
Analyzing the transition layer between two phases we estimate the surface tension coefficient and 
speculate on the possible existence of quark droplet. 
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Notwithstanding the well-known incompletion of quan- 
tum chromodynamics (QCD) and its strongly limited ca- 
pacity to perform the non-perturbative calculations both 
in the vacuum and at the finite temperature and bary- 
onic density, one of the major predictions of this theory, a 
possible existence of quark-gluon plasma, was so exciting 
and convincing that pushed forward very active research 
programme in experiments with relativistic heavy ions. 
In such a situation when the applications which should 
usually be based on the standard interrelations between 
the hadronic properties and the QCD Lagrangian param- 
eters are merely impossible, but the practical need in the 
quantitative estimates is dictated by the running experi- 
ments, we are forced to be very persistent and pragmatic 
in searching the effective Lagrangians. 

The Nambu-Jona-Lasinio (NJL) model and its numer- 
ous extensions are the most popular in this context be- 
cause they share some global symmetries of QCD and 
allow us to make surmountable some serious difficulties 
and uncertainties faced in the QCD calculations. It ap- 
pears especially appreciable and effective in studying the 
nature of nuclear matter and its (super) dense state be- 
ing treated as a model of QCD at large quark chemical 
potential. Nowadays the experiments in heavy-ion colli- 
sions to a considerable extent are driven by the results 
of phenomenological investigations of the properties of 
nucleon-nucleon force and the phase diagram of strongly 
interacting matter which relies on the corresponding es- 
timates of experimentally measurable quantities. 

Thus, exploring the QCD phase diagram with the ef- 
fective models is targeted from the theoretical view point 
by necessity to find out some kind of interpolation be- 
tween the physics as conceived by the lattice QCD sim- 
ulations (still unrealistic because of several reasons) and 
the physics outputs of phenomenological studies. The 



*Also at Institute of Theoretical and Experimental Physics, 
Moscow, RUSSIA 



vast activity (and progress) along this line [l| together 
with the recent experimental results from LHC and RHIC 
lead to the unexpected vision of still pending questions 
and, perhaps, their new interpretation. 

These and some related topics are discussed in this pa- 
per inspired by well known and fruitful idea about the 
specific role of surface degrees of freedom in the finite 
fermi-liquid systems d and, to a considerable extent, 
by our previous works 4^, [sl in which the quarks were 
treated as the quasi-particles of the model Hamiltonian, 
and the problem of filling up the Fermi sphere was studied 
in detail. In particular, within the Nambu-Jona-Lasinio 
(NJL) model [6] new solution branches of the equation 
for dynamical quark mass as a function of chemical po- 
tential (the details are shown below in Fig. [IJ have been 
found out. Besides, the existence and origin of the state 
filled up with quarks which is almost degenerate with the 
vacuum state both in the quasi-particle chemical poten- 
tial and in the ensemble pressure has been demonstrated. 
In general, the approach developed may be considered as 
another microscopical substantiation of the bag model in 
which the states filled up with quarks might be instru- 
mental 'construction material' for baryons. 

Our analysis here is performed within two approaches 
which are supplementary, in a sense, but, fortunately, 
lead to the identical results. One of these approaches 
based on the Bogolyubov transformation is especially in- 
formative for studying the process of filling the Fermi 
sphere up because the density of quark ensemble devel- 
ops a continuous dependence on the Fermi momentum in 
this case. It allows us to reveal an additional structure 
in the solution of gap equation for dynamical quark mass 
just in the proper interval of parameters characteristic 
for the phase transition and to trace its evolution. It 
results in the possibility for quark (fermionic) ensemble 
to be found in two aggregate states, a gas and a liquid, 
and the partial restoration of chiral condensate in a liq- 
uid phase (Section I). In order to make these conclusions 
easily perceptible we deal with the simplest version of the 
NJL model (with one flavor and one of the standard pa- 
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rameter sets). We try also to construct a description of 
transition layer between two phases and, in particular, to 
estimate the surface tension coefRcient (Section II) what 
is of obvious importance in context of discussing the pos- 
sible quark droplet formation (Section III). Some tech- 
nical moments of calculating the mean energy functional 
are picked out into Appendix. 



I. EXPLORING QUARK ENSEMBLE 

Now as an input for starting we remind the key el- 
ements of approach which has been developed 
The corresponding model Hamiltonian includes the inter- 
action term taken in the form of a product of two colour 
currents located in the spatial points x and y which are 
connected by a form-factor and its density reads as 



U 



-q{i^W + im)q-qe^^q / dyq't'^^q' {A^X) , (1) 



where q = q(x), q = ^(x), q' = q(y), q' = g(y) are the 
quark and anti-quark operators. 



7 (2^)3 (2^)1/^ Kp,s,c)u„,(p,s,c)e*P^- 
-h6+(p,s,c)w„i(p,s,c)e"*P''] , 



(2) 



p\ — — — TO^, i-is the colour index, a is the spinor 
index in the coordinate space, a^, a and h are the 
creation and annihilation operators of quarks and anti- 
quarks, a |0) = 0, 6 |0) — 0, |0) is the vacuum state of 
free Hamiltonian and to is a current quark mass. The 
summation over indices s and c is meant everywhere, 
the index s describes two spin polarizations of quark and 
the index c plays the similar role for a colour. As usual 
= A"/2 are the generators of SU{Nc) colour gauge 
group. The Hamiltonian density is considered in the Eu- 
clidean space and 7^ denote the Hermitian Dirac matri- 
ces, = 1,2,3,4. {AIj^A'^) stands for the form-factor 
of the following form 
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Y [^(x ~ y)'^M'> ~ Jp-i^ ~ y)] ' (3) 



where the second term is spanned by the relative distance 
vector and the gluon field primed denotes that in the 
spatial point y. The effective Hamiltonian density ([T]) 
results from averaging the ensemble of quarks influenced 
by intensive stochastic gluon field see Ref. 01 . For 
the sake of simplicity we neglect the contribution of the 
second term in ([3]) in what follows. The ground state of 
the system is searched as the Bogolyubov trial function 
composed by the quark-anti-quark pairs with opposite 
momenta and with vacuum quantum numbers, i.e. 



k) = r |o) , 

r = Hp,,, exp{.^[a+(p, s)fe+(-p, s) + a(p, s)b{~p, s)]}' 



(4) 



In this formula and below, in order to simplify the nota- 
tions, we refer to one compound index only which means 
both the spin and colour polarizations. The parameter 
</j(p) which describes the pairing strength is determined 
by the minimum of mean energy 



E=MH\a) 



(5) 



By introducing the 'dressing transformation' we define 
the creation and annihilation operators of quasi-particles 
'AS A ~ T a T^^, = T b^T^^ and for fermions 
T^^ — ■ Then the quark field operators are presented 
as 



9(x) 



9(x) 



dp 1 



(2^)3 (2|p4|)i/ 
-l-B+(p,s) y(p,s) e-^P''] , 
dp 1 



— [A{p,s)U{p,s)e^^-- 



— [A+(p,s) C/(p,s) e-'P^+ 



(2^)3|2|p4|)i/ 
+B(p,s)F(p,s) e'P-] , 

and the transformed spinors U and V are given by the 
following forms 



J7(p, s) ~ cos{ip) u{p, s) — sin((^) w(— p, s) 
y(p, s) = sin(<y9) u(-p, s) + cos{^p) v{p, s) 



(6) 



where C/(p, s) — U^{p,s) 74, V{p,s) — V^{p,s) 74 are 
the Dirac conjugated spinors. 

In Ref. fB*! the process of filling in the Fermi sphere 
with the quasi-particles of quarks was studied by con- 
structing the state of the Sletter determinant type 



\N) 



(7) 



\P\<Pf;S 



which possesses the minimal mean energy over the state 
|A^). The polarization indices run through all permissi- 
ble values here and the quark momenta are bounded by 
the limiting Fermi momentum Pp. The momenta and 
polarizations of states forming the quasi-particle gas are 
marked by the capital letters similar to above formula 
and the small letters are used in all other cases. 

As it is known the ensemble state at finite tempera- 
ture T is described by the equilibrium statistical opera- 
tor ^. Here we use the Bogolyubov-Hartree-Fock approx- 
imation in which the corresponding statistical operator 
is presented by the following form 



-P He, 



Tr {e 



-0 H. 



app 



} , 



(8) 



where an approximating effective Hamiltonian iJapp is 
quadratic in the creation and annihilation operators of 
quark and anti-quark quasi-particles A"*", A, B~^, B and 
is defined in the corresponding Fock space with the vac- 
uum state I cr) and /3 = T-i. There is no need to know 
the exact form of this operator henceforth because all 
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the quantities of our interest in the Bogolyubov-Hartree- 
Fock approximation are expressed by the correspond- 
ing averages n(P) = Tr{^A+(P; 5')A(P; S*)}, n{Q) = 
Tr{^B+(Q;r)B(Q;T)} which are obtained by solving 
the following variational problem. The statistical oper- 
ator ^ is defined by such a form in order to have the 
minimal value of mean energy of quark ensemble 

E = Tr{C H} 

at the fixed mean charge 

dp 



]i = Tr{e Qi} = V 2N, 



(27r)3 



[nip) - n{p)] , (9) 



where 



Qa = 



— I d:x. qi^iq 

dp —ip4 
(27r)3 \p,\ 



[A+{p)A{p) + B{p)B+{p)] 



for the diagonal component (which is a point of our in- 
terest here) and at the fixed mean entropy (S — — In^) 



S 



TrUln^^ 
dp 



-V2N, 



(27r)3 
+n(ji) lnh{p) + (1 



[n(jj) In n{p) - 
— n{p)) ln(l 



(1 - n{p)) ln(l 
- nip))] ■ 



(10) 

n{p)) + 



The mean charge ^ is calculated here up to the unessen- 
tial (infinite) constant coming from permuting the op- 
erators BB^ in the charge operator Q4. It is appro- 
priate here to remind that the mean charge should be 
treated in some statistical sense because it characterizes 
quark ensemble density and has no colour indices. The 
mean energy density per one quark degree of freedom 
w = £/{2Nc), £ = E/V where E \s & total energy of 
ensemble is calculated (the details of derivation can be 
found in Appendix) to get the following form 



dp 

(2^ 



|P4| + 



-^^|P4| cos9[n+n 

[ZTT j'' 



-n-aj 



dp 



sin (6 - 



\n+n 



dq 

(2^ 



■ sm I 



9'-e'J[n'+n'-l]I, (11) 



r 



where = 2(p, 6' = e{q), n' = n{q), I = /(p + q) and the 
angle 0rn{p) is determined by sin^m — m/\p4\. We are 
interested in minimizing the following functional 



n^E-^Q4-TS 



(12) 



where fj, and T are the Lagrange factors for the chemical 
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potential and temperature respectively. The approximat- 
ing Hamiltonian H^pp is constructed simply by using the 
information on E — fi Q4 of presented functional (see, 
also below). For the specific contribution per one quark 
degree of freedom / = F/{2Nc), F ~ il./V we receive 



/ = 



dp 

(2^ 
rfq 



[Ipil cosO (n + n — 1) ~ fi [n — n)] + 



dp 

J2nf- 



G 



dp 

(2^ 



sin [9 — 9m) [n + n — 1) x 



(27r)3 



9'^) {n' + n'-l)I + T 



dp 

(2^ 



(13) 



[n In n -|- (1 — n) ln(l — n) -f n In n + (1 — n) ln(l — n)] 



The optimal values of parameters are determined by solv- 
ing the following equation system {df / d9 = 0, df / dn ~ 0, 
df/dn = 0) 

|P4| sin6'- Mcos(6'- 6'„) = , (14) 
|P4| cos 61 - /I -I- M sin (61 - 9^) - T In (n^^ - l) = , 
\pi\ cos9 + n + M sin (6* - 9„,) - T In {n.-^ - l) = , 

where we denoted the induced quark mass as 
dq 



M(p) = 2G, 



(2^)= 



-(l-n'-rl') sin {9' - 9'J /(p+q). (15) 



Turning to the presentation of obtained results in the 
form customary for the mean field approximation we in- 
troduce a dynamical quark mass Mq parametrized as 



sm 



^ |P4| = (p^ + M|(p))V2, (ig) 



and ascertain the interrelation between induced and dy- 
namical quark masses. From the first equation of system 
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FIG. 1: The residual A for equation fTS} is presented as a 
function of dynamical quark mass Af, (MeV) at zero value of 
temperature and the following values of chemical potential fj, 
(MeV) - 335 (the lowest curve), 340, 350, 360 (the top curve). 



we fix the pairing angle 



smf 



p M 

\P4m 



and making use the identity 



(b4 



M mf + M^p"^ = [p2 + (M - m)2] \p^\^ , (17) 



find out that 



cos 9 = ± 



\Pi 



m M 



\p4\p4 



For the sake of clarity we choose the upper sign 'pins'. 
Then, as an analysis of the NJL model teaches, the 
branch of equation solution for negative dynamical quark 
mass is the most stable one. Let us remember here we 
are dealing with the Euclidean metrics (though it is not 
a principal point) and a quark mass appears in the corre- 
sponding expressions as an imaginary quantity. Now sub- 
stituting the calculated expressions for the pairing angle 
into the trigonometrical factor 



sm ( 



sin 9 



P 



\P4 



— COSf 



m 

R 



and performing some algebraic transformations we come 
to the relation 



(18) 



M,(p) = M(p) - m . 

In particular, the equation for dynamical quark mass (|15| 
is getting the form characteristic for the mean field ap- 
proximation 



M = 2G 



(1 



n')^/(p + q). (19) 



The second and third equations of system ([T^ allow 
us to find the following expressions 



(20) 



FIG. 2: The chemical potential fi (MeV) is plotted as a func- 
tion of charge density Q4 — Q4/{3V) (in the units of ch/fm^). 
The factor 3 relates the densities of quark and baryon matters. 
The top curve corresponds to the situation of zero tempera- 
ture. The curves following down correspond to the temper- 
ature values r = 10 MeV, ... , T = 50 MeV with spacing 
T = 10 MeV. 



and, hence, the thermodynamic properties of our system, 
in particular, the pressure of quark ensemble 

dE 
dV 



P= — ^ . 



By definition we should calculate this derivative at con- 
stant mean entropy, dS/dV = 0. This condition makes 
possible, for example, to calculate the derivative dfi/dV, 
but the mean charge Q4 should not also change. In or- 
der to maintain it valid we introduce two independent 
chemical potentials — for quarks /x and for anti-quarks p, 
(following Eq. ((20| with the opposite signs) . It leads also 
to the change /i — ;> /i in definition of n in Eq. (refnew- 
den). This kind of description apparently allows us to 
treat even some non-equilibrium states of quark ensem- 
ble (but with losing a covariance similar to the situation 
which takes place in electrodynamics while one deals with 
electron-positron gas). Here we are interested in the un- 
affected balanced situation of /i = fj,. Then the corre- 
sponding derivative of specific energy dw/dV might be 
presented as 



dw 
dV 



dp f dn dfi dn dp. 



(27r)3 \dfidV dpdV 



[\P4\ C0S( 



-2Gsin{ 



sin {9' - 9[^) {n' + n' - 1)1 



Now representing the trigonometric factors via dynami- 
cal quark mass and drawing Eq. (|15p we obtain for the 
ensemble pressure 



P = 



E 

'v 



dp f dn dp dn d/j, 



(27r)3 \dfidV ' dfldV 
The requirement of mean charge conservation 



\P,\. (21) 



dV V 



dp / dn dp dn dfi . 
(27r)3 [jildV^dfidV " ' 



(22) 
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FIG. 3: The ensemble pressure P (MeV/fm"^) is shown as a 
function of charge density Q4 at temperatures T — MeV, ... 
, T = 50 MeV with spacing T = 10 MeV. The lowest curve 
corresponds to zero temperature. The dashed curve shows 
the boundary of phase transition liquid-gas, see the text. 
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FIG. 4: The fragments of isotherms in Fig. (2] |3l see text. 
Chemical potential (MeV) is plotted as a function of pres- 
sure P MeV/fm^. The top curve corresponds to the zero 
isotherm and following down with spacing 10 MeV till the 
isotherm 50 MeV (the lowest curve). 



provides us with an equation which interrelates the 
derivatives dn/dV and djl/dV. Apparently, the regu- 
larized expressions for mean charge of quarks and anti- 
quarks are meant ([9]) here. Dealing in a similar way with 
the requirement of mean entropy conservation, dS / dV — 
0, we receive another equation as 



dp dn 



In- 



dfi 



(27r)3 d/i 1 
S 

2Nc ' 



dV 



dp dn 



In- 



d^ 



(27r)3 dp. 1 



^dV 



(23) 



Substituting here Tln(n~-'^ — 1) = — /i -I- IP4I and 
rin(n~^ — 1) — /H- IP4I we have after simple calculations 
with taking into account that 



dp / dn dfi dn dfi 



(27r)3 \d^idV dfidV 



\Pa 



ST 



Eventually it leads to the following expression for the 
pressure 



P 



E 
V 



S T 



V 



(24) 



(of course, the thermodynamic potential is 57 = —P V). 
At small temperatures the anti-quark contribution is neg- 
ligible, and thermodynamic description can be grounded 
on utilizing one chemical potential /i only. If the anti- 
quark contribution is getting intrinsic the thermody- 
namic picture becomes more complicated due to the ne- 
cessity to obey the condition p. = fJ. which comes to the 
play. In particular, at zero temperature we might con- 
sider the anti-quark contribution absent and obtain 



P = -£ 



A* Pq 



1 /2 

where ^ = [P^ + M^{Pf)) , Pf is the Fermi momen- 
tum and Pq — N/V is the quark ensemble density. 



For lucidity of our view point we consider mainly the 
NJL model [6] in this paper, i.e. the correlation function 
(the form- factor in Eq. ([3])) behaves as the ^- function in 
coordinate space. It is well known fact that in order to 
have an intelligent result in this model one needs to use 
a regularization cutting of the momentum integration in 
Eq. (|13p . We adjust the standard set of parameters [7| 
here with |p| < A, A = 631 MeV, m = 5.5 MeV and 
GAV(27r2) = 1.3. This set at n = 0, n = 0, T = 
gives for the dynamical quark mass Mq = 335 MeV. In 
particular, it may be shown the following representation 
of ensemble energy is valid at the extremals of functional 

m 



E = E,ac + 2N,V 



dp 

(2^ 



IP4I (n + n) 



'^'P .{\p^\-\P^\) + 2N,V^, 



(25) 



(27r)3^'""' ' AG' 

It is easy to understand this expression with the vacuum 
contribution subtracted looks like the energy of a gas of 
relativistic particles and antiparticles with the mass Mq 
and coincides identically with that calculated in the mean 
field approximation. 

Let us summarize the results of this exercise. So, we 
determine the density of quark n and anti-quark n quasi- 
particles at given parameters /i and T from the second 
and third equations of system (fHl) . From the first equa- 
tion we receive the angle of quark and anti-quark pairing 
as a function of dynamical quark mass Mq which is han- 
dled as a parameter. Then at small temperatures, below 
50 MeV, and value of chemical potentials of dynamical 
quark mass order, p, ~ Mq, there are several branches of 
solutions of the gap equation. Fig. [1] displays the differ- 
ence of right and left sides of Eq. (|15p which is denoted 
by A at zero temperature and several values of chemical 
potential ^ (MeV) = 335 (the lowest curve), 340, 350, 
360 (the top curve) as a function of parameter Mq. The 
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FIG. 5: The dynamical quark mass \Mq\ (MeV) is shown as 
a function of chemical potential ^ (MeV) at the temperatures 
T = MeV, ... , T = 100 MeV with spacing T = 10 MeV. 
The most right curve corresponds to zero temperature. 



zeros of function A(Mg) correspond to the equilibrium 
values of dynamical quark mass. 

The evolution of chemical potential as a function 
of charge density Q4 = Qi/{'iV) (in the units of 
charge/ /m'^) with the temperature increasing is depicted 
in Fig. [2] (factor 3 relates the quark and baryon matter 
densities). The top curve corresponds to the zero temper- 
ature. The other curves following down have been calcu- 
lated for the temperatures T = 10 MeV, ... , T = 50 MeV 
with spacing T = 10 MeV. As it was found in Ref. [ll 
the chemical potential at zero temperature is increasing 
first with the charge density increasing, reaches its maxi- 
mal value, then decreases and at the densities of order of 
normal nuclear matter density (l9|. Pq ~ 0.16//m"^, be- 
comes almost equal its vacuum value. Such a behaviour 
of chemical potential results from the fast decrease of 
dynamical quark mass with the Fermi momentum in- 
creasing. It is clear from Fig. [5] the charge density is 
still multivalued function of chemical potential at the 
temperature slightly below 50 MeV. The Fig. [3] shows 
the ensemble pressure P (MeV/fm'^) as the function of 
charge density Q4 at several values of temperature. The 
lowest curve corresponds to the zero temperature. The 
other curves following up correspond to the temperatures 
T = 10 MeV, ... , T = 50 MeV (the top curve) with spac- 
ing T = 10 MeV. It is curious to remember now that in 
Ref. [5I the vacuum pressure estimate for the NJL model 
was received as 40—50 MeV/fm^^ which is entirely com- 
patible with the results of conventional bag model. Be- 
sides, some hints at an instability presence (rooted in the 
anomalous behavior of pressure dP/dn < 0, see also Q, 
[^) in some interval of the Fermi momentum has been 
found out. 

Fig. [3] shows the fragments of isotherms of Fig. [51 [3] 
but in the different coordinates (chemical potential — en- 
semble pressure) . The top curve is calculated at the zero 
temperature, the other isotherms following down corre- 
spond to the temperatures increasing with spacing 10 
MeV. The lowest curve is calculated at the temperature 
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T 

FIG. 6: The dynamical quark mass \Mq\ (MeV) as a function 
of temperature at the small value of charge density Q4. 



50 MeV. This plot obviously demonstrates that there are 
the states on the isotherm which are thermodynamically 
equilibrated and have an equal pressure and chemical 
potential (see the characteristic Van der Waals triangle 
with the crossing curves). The equilibrium points cal- 
culated are shown in Fig. [3] by the dashed curve. The 
points of dashed curve crossing with an isotherm pin- 
point the boundary of gas — liquid phase transition. The 
corresponding straight line P = const which obeys the 
Maxwell rule separates the non-equilibrium and unsta- 
ble fragments of isotherm and describes a mixed phase. 
Then the critical temperature for the parameter which 
we are using in this paper becomes Tc ^ 45.7 MeV with 
the critical charge density as Q4 ~ 0.12 ch/fm'^. Usu- 
ally the thermodynamic description is grounded on the 
mean energy functional which is the homogeneous func- 
tion of particle number like E ^ N /{S/N, V/N) (with- 
out vacuum contribution). It is clear such a description 
requires the corresponding subtractions to be introduced, 
however, this operation does not change the final results 
considerably. Now the vague arguments of Refs. Q that 
the states filled up with quarks and separated from the 
instability region look like a 'natural construction mate- 
rial' to form the baryons are getting much more clarity 
and give a hope to understand the existing fact of equi- 
librium between vacuum and octet of stable (in strong 
interaction) baryons[20]. 

The dynamical quark mass \Mq\ (MeV) as a function of 
chemical potential (MeV) is presented for the tempera- 
tures T = MeV, ... , T = 100 MeV with spacing T = 10 
MeV in Fig. [S] The most right-hand curve corresponds to 
the zero temperature. At small temperatures, below 50 
MeV, the dynamical quark mass is the multivalued func- 
tion of chemical potential. Fig. [6] shows the dynamical 
quark mass as a function of temperature at small values 
of charge density Q4 ~ 0. Such a behaviour allows us 
to conclude that the quasi-particle size is getting larger 
with temperature increasing. It becomes clear if we re- 
member that the momentum corresponding the maximal 
attraction between quark and anti-quark pg (according 
to Ref. [4]) is defined by dsinO/dp = 0. In particular. 
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this parameter in the NJL model equals to 

pe = (|M,| m)i/2 . (26) 

but its inverse magnitude defines the characteristic (ef- 
fective) size of quasi-particle re — pg^ . 

If one is going to define the quark chemical potential as 
an energy necessary to add (to remove) one quasi-particle 
(as it was shown in [5] at zero temperature), /i = dE/dN, 
then in vacuum (i.e. at quark density pg going to zero) 
quark chemical potential magnitude coincides with the 
quark dynamical mass. It results in the phase diagram 
displayed at this value of chemical potential although, 
in principle, this value could be smaller than dynamical 
quark mass as it has been considered in the pioneering 
paper [10]. If one takes, for example, chemical potential 
value equal to zero it leads to the conventional picture 
but, obviously, such a configuration does not correspond 
to the real process of filling up the Fermi sphere with 
quarks. 

Apparently, our study of the quark ensemble thermo- 
dynamics produces quite reasonable arguments to pro- 
pound the hypothesis that the phase transition of chi- 
ral symmetry (partial) restoration has been already real- 
ized as the mixed phase of physical vacuum and baryonic 
matter (2H. However, it is clear our quantitative esti- 
mates should not be taken as appropriate for compar- 
ing with, for example, the critical temperature of nuclear 
matter phase transition which has been experimentally 
measured and is equal 15-20 MeV. Besides, the gas com- 
ponent (at T = 0) has nonzero density (0.01 of the nor- 
mal nuclear density) but in reality this branch should 
correspond to the physical vacuum, i.e. zero baryonic 
density [22'|. In principle, an idea of global equilibrium 
of gas and liquid phases pursued us to put the adequate 
boundary conditions down at describing the transitional 
layer existing between the vacuum and the filled-up state 
and to calculate the surface tension effects. It looks plau- 
sible that the changes taking place in this layer could 
ascertain all ensemble processes similar to the theory of 
Fermi-liquids. 

II. TRANSITION LAYER BETWEEN GAS AND 
LIQUID 

This concept advanced would obtain the substantial 
confirmation if we are able to demonstrate an existence 
of this transition layer at which the ensemble transfor- 
mation from one aggregate state to another takes place. 
As it was argued above the indicative characteristic to 
explore a homogeneous phase (at finite temperature) is 
the mean charge (density) of ensemble. It was demon- 
strated the other characteristics, for example, a chiral 
condensate, a dynamical quark mass, etc. could be recon- 
structed as well. So, here we are analyzing the transition 
layer at zero temperature. 

If one assumes the parameters of gas phase are approx- 
imately the same as those at zero charge pg = 0, i.e. as 



in vacuum (it means ignoring the negligible distinctions 
in the pressure, chemical potential and quark conden- 
sate), the dynamical quark mass develops the maximal 
value, and for the parameter choice of the NJL model 
it is M = 335 MeV. Then from the Van der Waals 
diagram one may conclude that the liquid phase be- 
ing in equilibrium with the gas phase develops the den- 
sity pi = 3 X 0.185 ch/fm'^ (by some reason which be- 
comes clear below we correct it in favour of the value 
pi = 3 X 0.157 ch/fm"^). The detached factor 3 here links 
again the magnitudes of quark and baryon densities. The 

quark mass is approximately M~ 70 MeV in this phase 
(and we are dealing with the simple one-dimensional pic- 
ture hereafter). 

The precursor experience teaches that an adequate de- 
scription of heterogeneous states can be reached with 
the mean field approximation [l^ . In our particu- 
lar case it means making use the corresponding ef- 
fective quark-meson Lagrangian fisj (the functional of 
Ginzburg-Landau type) 

C = -q (d + M) q - i {d^a)' - U{a) - 

(27) 

1 771 

F^^Ff,^ - V^Vf, - ga qq (T + ig^ g 7m 9 ' 

where 

F^. = - d,V^ , Uia) = !!^ a' + ^a' + ^a\ 

and a is the scalar field, is the field of vector mesons, 
TOct, ruy are the masses of scalar and vector mesons and 
go-, gv are the coupling constants of quark-meson inter- 
action. The U (cr) potential includes the nonlinear terms 
of sigma field interactions up to the fourth order. For 
the sake of simplicity we do not include the contribution 
coming from the pseudoscalar and axial- vector mesons. 

The meson component of such a Lagrangian should be 
selfconsistently treated by considering the corresponding 
quark loops. (In terms of a relativistic extension of the 
Landau theory of Fermi-liquid the density fluctuations 
(meson field collective modes) are nothing more than the 
zero-sound as was shown in Ref. [9]). Here we do not see 
any reason to go beyond the well elaborated and reliable 
one loop approximation (P7)l [Tsl, although recently the 
considerable progress has been reached (as we mention at 
the beginning of this paper) in scrutinizing the nonhomo- 
geneous quark condensates by application of the powerful 
methods of exact integration [1J|. Here we believe it is 
more practical to adjust phenomenologically the effective 
Lagrangian parameters basing on the transparent physi- 
cal picture. It is easy to see that handling (|27l) in one loop 
approximation we come, in actual fact, to the Walecka 
model [15] but adopted for the quarks. In what follows 
we are working with the designations of that model and 
do hope it does not lead to the misunderstandings. 

In the context of our paper we propose to interpret 
Eq. (1271) in the following way. Each phase might be con- 
sidered, in a sense, with regard to another phase as an 
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excited state which requires the additional (apart from a 
charge density) set of parameters (for example, the meson 
fields) for its complete description, and those are char- 
acterizing the measure of deviation from the equilibrium 
state. Then the crucial question becomes whether it is 
possible to adjust the parameters of effective Lagrangian 
(I27p to obtain the solutions in which the quark field inter- 
polates between the quasi-particles in the gas (vacuum) 
phase and the quasi-particles of the filled-up states. For 
all that the density of the filled-up state ensemble should 
asymptotically approach the equilibrium value of pi and 
should turn to the zero value in the gas phase (vacuum) . 

The scale inherent in this problem may be assigned 
with one of the mass referred in the Lagrangian ((27|). In 
particular, we bear in mind the dynamical quark mass in 
the vacuum M . Besides, there are another four indepen- 
dent parameters in the problem and in order to compare 
them with the results of studying a nuclear matter we 
employ the form characteristic for the (nuclear) Walecka 
model 



M 



M 



9l M 



9. 



Taking the parameterization of the potential U{a) as 
ba = 1.5 ml {ga/M), Ca = 0.5 ml (gcr/M)"^ we come 
to the sigma model but the choice 6 = 0, c = results in 
the Walecka model. As to standard nuclear matter appli- 
cation the parameters b and c demonstrate vital model 
dependent character and are quite different from the pa- 
rameter values of sigma model. Truly, in that case their 
values are also regulated by additional requirement of an 
accurate description of the saturation property. On the 
other hand, for the quark Lagrangian (|27p we could intu- 
itively anticipate some resemblance with the sigma model 
and, hence, could introduce two dimensionless parame- 
ters 77 and ^ in the form oi b = rj ba, c = C,'^ Ca which 
characterize some fluctuations of the effective potential. 
Then the scalar field potential is presented as follows 



m^ 0^ ( 



4 — rya-fCV) ^2 



The meson and quark fields are defined by the following 
system of the stationary equations 

A cr — TO^ G = b ^ c ^ ga Ps , 

~mlV ^-g, p , (28) 

(V+ M) q={E-g,Y)q, 

where M— A/ + g„G is the running value of dynamical 
quark mass, E stands for the quark energy and V = 
—iVi- The density matrix describing the quark ensemble 
at T = has the form 



e(^) 



Pf 



dp 

(2^ 



qp{x) qp{x) , 



in which p is the quasi-particle momentum and the Fermi 
momentum Pp is defined by the corresponding chemical 



potential. The densities ps and p at the right hand sides 
of Eq. ([^5]) are by definition 



p,(x) -TrU(x),l} , p{x)=Tr{ax),'r4 . 

Here we confine ourselves to the Thomas-Fermi ap- 
proximation while describing the quark ensemble. Then 
the densities which we are interested in are given with 
some local Fermi momentum Pf{x) as 



P = l 



dp 

(27r)3 

dp 



7 
67r2 



F J 



(29) 



_ M 

(27r)3 E 



47r^ 



2\l/2 



In 



(1 + A^) 



1/2 



(l + A2)i/2-f 



where 7 is the quark factor which for one flavour is 7 = 
2Nc {Nc is the number of colours), E = (p^-f M 

and X ~M /Pp- By definition the ensemble chemical 
potential does not change and it leads to the situation in 
which the local value of Fermi momentum is defined by 
the running value of dynamical quark mass and vector 
field as 



p = M ^ g^V + 



M 



1/2 



(30) 



Now we should tune the Lagrangian parameters ([?f)) . 
For asymptotically large distances (in the homogeneous 
phase) we may neglect the gradients of scalar and vector 
fields and the equation for scalar field in the system (pS)) 
leads to the first equation which bounds the parameters 
Cs J C'w , b, c 



M\M " M) 



-b A.]{M - Mf+c[M - Mf 



'Ps. (31) 



The asymptotic vector field is given by the ensemble den- 
sity V ~ Cy p/ {gvAP). The second equation results from 
the relation (|30p for the chemical potential and gives 



M = 



C,, 



M2 



Pf+ M 



1/2 



(32) 



Extracting the liquid density from ([^H) we obtain the 
Fermi momentum (Pp = 346 MeV) and applying the 
identities (PT|) . (|32p we have for the particular case 6 = 0, 
c = that = 25.3, = -0.471, i.e. the vector com- 
ponent is small (comparing to Cf) and has negative 
value which is unacceptable. Apparently, it looks neces- 
sary to neglect the contribution coming from the vector 

field or to diminish the dynamical quark mass M up to 
the value which retains the identity (I32p valid with posi- 
tive or equals to zero. In the gas phase the dynamical 
quark mass can also be corrected to the value larger than 
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FIG. 7: The domain of the 77, c ((" = c 7;)-plane in which an 
increase of specific energy occurs, see the text. The dots rep- 
resent a stable kink. The star shows the position of canonical 
(chiral) kink, see the text. 



FIG. 8: The stable kink solutions with c — 1.1, the solid line 
corresponds to rj ~ 0.977 {rria ~ 468 MeV) and the dashed 
line corresponds to « 1.813 {rricr ~ 690 MeV), x is given in 
the units of fm and a is given in MeV. 



the vacuum value. It is clear that in the situation of the 
hquid with the density p/ = 3 x 0.185 ch/fm"^ the dynam- 
ical quark mass should coincide (or exceed) M = 346 
MeV in the gas phase. However, here we correct the liq- 
uid density (as it was argued above) to decrease its value 
up to /O/ = 3 X 0.157 ch/fm'^ which is quite acceptable in 
the capacity of normal nuclear matter density. In fact, 
this possibility can be simply justified by another choice 

of the NJL model parameters. Thus, we obtain at M= 70 
MaB and 6 = 0, c = that = 28.4, = 0.015, i.e. 
we have a small but positive value for the vector field 
coefficient. At the same time, being targeted here to es- 
timate the surface tension effects only we do not strive 
for the precise fit of parameters. In the Walecka model 
these coefficients are = 266.9, = 145.7, (6 = 0, 
c — 0). Moreover, there is another parameter set with 
Cg = 64., ss |l3| but it is rooted in an essential non- 
linearity of the sigma-field due to the nontrivial values of 
the coefficients b and c. The option (formally unstable) 
with negative c (6) has been also discussed. 

The coupling constant of scalar field is fixed by the 
standard (for the NJL model) relation between the quark 
mass and the 7r-meson decay constant g^r = M / fj^ (we 
put /tt = 100 MeV) although there is no any objection to 
treat this coupling constant as an independent parame- 
ter. As a result of all agreements done we have for the a- 
meson mass ma- — M/Cg- In principle, we could even 
fix the CT-meson mass and coupling constant 170. but all 
relations above mentioned lead eventually to quite suit- 
able values of the cr-meson mass as will be demonstrated 
below. 

The vector field plays, as we see, a secondary role be- 
cause of the small magnitude of constant Cy . Then tak- 
ing the vector meson mass as ~ 740 MeV (slightly 
smaller than the mass of w-meson because of simple tech- 
nical reason only) we calculate the coupling constant of 
vector field from the relation similar to the scalar field 
TO„ = M/Cy. Amazingly, its value comes about 



steadily small in comparing to the value characteristic 
for the NJL model g^ — \/&ga- However, at the values 
of constant which we are interested in it is very diffi- 
cult to maintain the reasonable balance and for the sake 
of clarity we prefer to choose the massive vector field. 
Actually, it is unessential because we need this param- 
eter (as we remember) only to estimate the vector field 
strength. 

The key point of our interest here is the surface tension 
coefficient [16] which can be defined as 



= At: ri 



dx 



£{x)~— p{x) 
Pi 



(33) 



The parameter will be discussed in the next section at 
considering the features of quark liquid droplet, and for 
the present we would like to notice only that for the pa- 
rameters considered its magnitude for Nf = 1 is around 
To = 0.79 fm. Recalling the factor 3^/'^ which connects 
the baryon and quark numbers one can find the magni- 
tude (ro — 3^/'^0.79 « 1.14 fm) is in full agreement with 
the magnitude standard for the nuclear matter calcula- 
tions (in the Walecka model) = 1.1 — 1.3 fm. 

In order to proceed we calculate £{x) in the Thomas- 
Fermi approximation as 



£{x) 



7 



Pf{x) 



(27r)3 



M (x) 



1/2 



+ ^ 9v P{x) V{x) 9a Ps[x) (j{x) 



And to give some idea for the 'setup' prepared we present 
here the characteristic parameter values for some fixed h 
andcwith/3; = 3x0.157ch/fm'^. In the liquid phase they 

are M= 70 MeV {Pp = 327 MeV) and ei 310.5 MeV 
(index I stands for a liquid phase and e(yx) ~ £{x)/ p{x) 
defines the density of specific energy). Both equations 
(Pljl and ([5^ are obeyed by this state. There exist the 

solution with larger value of quark mass M= 306 MeV, 
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{Pp = 135 MeV) (we have faced the similar situation 
in the first section deahng with the gas of quark quasi- 
particles) and e = 338 MeV ~ eg {eg is the specific energy 
in the gas phase) which satisfies both equations as well. 
The specific energy of this solution occurs to be larger 
than specific energy of the previous solution. It is worth- 
while to mention the existence of intermediate state cor- 
responding to the saturation point with the mass M — 95 
MeV, (Pp = 291 MeV) and e = 306 MeV. Obviously, 
it is the most favorable state with the smallest value of 
specific energy (and with the zero pressure of quark en- 
semble) , and the system can reach this state only in the 
presence of significant vector field. This state (already 
discussed in the first section) corresponds to the minimal 
value of chemical potential (T = 0) and can be reached 
at the densities typical for the normal nuclear matter. 
However, Eq. ([32l) is not valid for this state. 

Two another parameters 77, C are fixed by looking 
through all the configurations in which the solution of 
equation system p8|) with stable kink of the scalar field 
does exist and describes the transition of the gas phase 
quarks to the liquid phase. First, it is reasonable to scan 
the 7?, c (C = c ?7)-plane, in order to identify the domain 
in which the increase of specific energy £ ~ £1 pj pi < 
is revealed at running through all possible states which 
provide the necessary transition (without taking into ac- 
count the field gradients). In practice one need to follow 
a simple heuristic rule. The state with Pp ~ 1 MeV (i.e. 
e and the corresponding p). The state of characteristic 
liquid energy £1 (together with pi) should be compared 
at scanning the Lagrangian parameters 77 and c. Just this 
domain where they are commensurable could provide us 
with the solutions in which we are interested and FiglT] 
shows its boundary. The curve could be continued be- 
yond the value r] — 2.5 but the values of corresponding 
parameter 77 are unrealistic and not shown in the plot. 

We calculate the solution of equation system p8)) nu- 
merically by the Runge-Kutta method with the initial 
conditions (t{L) « 0, (t'{L) ^ imposed at the large 
distance L ^ t where t is a characteristic thickness of 
transition layer (about 2 fm). Such a simple algorithm 
occurs quite suitable if the vector field contribution is 
considered as a small correction (what just takes place in 
the situation under consideration) and is presented as 

V{x)^-^ r dz e-^-^^-'^^ p{z) , 
2to« J-l 

where the charge (density) p is directly defined by the 
scalar field. We considered the solutions including the 
contribution of the vector field as well and the corre- 
sponding results confirm the estimates obtained. 

Rather simple analysis shows the interesting solutions 
are located along the boundary of discussed domain. 
Some of those are depicted in Fig. [7] as the dots. Fig. 
[5] shows the stable kinks of cr-field with the parame- 
ter c — 1.1 for two existing solutions with rj w 0.977 
(to^ « 468 MeV) (sohd line) and 7/ « 1.813 {m„ « 690 
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C 

FIG. 9: The surface tension coefficient Us in MeV as a function 
of parameter c {C, = c rf) for the curve of stable kinks (with 
7?< 1.2). 



MeV) (dashed line). For the sake of clarity we consider 
the gas (vacuum) phase is on the right. Then the asymp- 
totic value of (T-field on the left hand side (cr sa 80 MeV) 

corresponds to M — 70 MeV. The thickness of transition 
layer for the solution with 7/ ~ 0.977 is i « 2 fm while for 
the second solution i sa 1 fm. 

Characterizing the whole spectrum of the solutions ob- 
tained we should mention that there exist another more 
rigid (chiral) kinks which correspond to the transition 
into the state with the dynamical quark mass changing 
its sign, i.e. M — 7> —M. In particular, the kink with the 
canonical parameter values 77 = 1, c = 1 is clearly seen 
(marked by the star in Fig. [7| and its surface tension 
coefficient is about 27n7r {m.^ is the 7r-meson mass). The 
most populated class of solutions consists of those hav- 
ing the meta-stable character. The system comes back 
to the starting point (after an evolution) pretty rapidly, 
and usually the tr-field does not evolve in such an extent 
to reach the asymptotic value (which corresponds to the 

dynamical quark mass in the liquid phase A/= 70 MeV). 
Switching on the vector field changes the solutions in- 
significantly (for our situation with small C^, it does not 
exceed 2 MeV in the maximum). 

The surface tension coefficient Ug in MeV for the curve 
of stable kinks with parameter 77 < 1.2 as the function of 
parameter c (^ = c 77) is depicted in Fig. 1^1 The a-meson 
mass at c « is 777^ ~ 420 MeV and changes smoothly 
up to the value 777o- w 500 MeV at c 1.16 (the maximal 
value of the coefficient c beyond which the stable kink so- 
lutions are not observed). In particular, 7710- ~ 450 MeV 
at c = 1. Two kink solutions with c ~ 1.1 for rj 0.977 
and for rj sa 1.813 (shown in Fig. [51 and the second one 
is not shown in Fig. [5]) have the tension coefficient values 
Us ~ 35 MeV and Us ~ 65 MeV, correspondingly. The 
maximal value of tension coefficient for the normal nu- 
clear matter does not exceed Us = 50 MeV. The nuclear 
Walecka model claims the value Tig w 19 MeV [16] as ac- 
ceptable and calculable. The reason to have this higher 
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FIG. 10: CT-field (MeV) as a function of the distance r (fm) 
for several solutions of the equation system (|28|l which are 
characterized by the net quark number Nq written to the left 
of each curve. 



FIG. 11: Distribution of the quark density p (ch/im^) for the 
corresponding solutions presented in Fig. 1111 



value of surface tension coefEcient for quarks is rooted 
in the different magnitudes of the mass deficit. Indeed, 

for nuclear matter it does not exceed M~ 0.5M albeit 
more realistic values are considered around M~ 0.7 M 
and for the quark ensemble the mass deficit amounts to 

M ~ 0.3M. We are also able to estimate the compres- 
sion coefficient of quark matter K which occurs signifi- 
cantly larger than the nuclear one. Actually, we see quite 
smooth analogy between the results of Section II and the 
results of bag soliton model [l3|- The thermodynamic 
treatment developed in the present paper allows us to for- 
mulate the adequate boundary conditions for the bag in 
physical vacuum and to diminish considerably the uncer- 
tainties in searching the true soliton Lagrangian. We be- 
lieve it was also shown here that to single out one soliton 
solution among others (including even those obtained by 
the exact integration method [ij])) which describes the 
transitional layer between two media, is not easy prob- 
lem if the boundary conditions above formulated are not 
properly imposed. 



III. DROPLET OF QUARK LIQUID 

The results of two previous sections have led us to put 
the challenging question about the creation and proper- 
ties of finite quark systems or the droplets of quark liquid 
which are in equilibrium with the vacuum state. Thus, 
as a droplet we imply the spherically-symmetric solution 
of the equation system (pS)) for a{r) and V{r) with the 
obvious boundary conditions cr'(O) = and V'{0) = in 
the origin (the primed variables denote the first deriva- 
tives in r) and rapidly decreasing at the large distances 
(7 —7- 0, F —7- 0, when r oo. 

A quantitative analysis of similar nuclear physics mod- 
els which includes the detailed tuning of parameters is 
usually based on the comprehensive fitting of available 
experimental data. This way is obviously irrelevant in 



studying the quark liquid droplets. This global difficulty 
dictates a specific tactics of analyzing. We propose to 
start, first of all, with selecting the parameters which 
could be worthwhile to play a role of physical observ- 
ables. Naturally, the total baryon number which phe- 
nomenologically (via factor 3) related to the number of 
valence quark in an ensemble is a reasonable candidate 
for this role. Besides, the density of quark ensemble p{r)^ 
the mean size of droplet Rq and the thickness of surface 
layer t look suitable for such an analysis. 

It is argued above that the vector field contribution is 
negligible because of the smallness of the coefficient Cy 
comparing to the Cs magnitude, and we follow this con- 
clusion (or assumption) albeit understand it is scarcely 
justified in the context of finite quark system. Thus, we 
will put gy — 0, V — in what follows and it will simplify 
all the calculations enormously. 

Fig. [TOl shows the set of solutions (cr-field in MeV) of 
the system (|28|) at Nf = 1 and Fig. [TT] presents the cor- 
responding distributions of ensemble density p (ch/fm'^). 
The parameters Cs , Cu , b and c are derived by the same 
algorithm as in the previous section, i.e. the chemical 
potential of quark ensemble M = 335 MsB (and a 0) 
is fixed at the spatial infinity. The fiUed-up states (liq- 

uid) are characterized by the parameters M— 70 MeV, 
Po = P/ = 3 X 0.157 ch/fm'^. The cr-meson mass and 
the coupling constant g^. are derived at fixed coefficients 
T] and ^, and they just define the behaviour of solutions 
(7(r), p{r), etc. The magnitudes of functions cr(r) and 
p{r) at origin are not strongly correlated with the values 
characteristic for the fiUed-up states and are practically 
determined by solving the boundary value problem for 
system (pS)) . In particular, the solutions presented in Fig. 
1101 have been received with the running coefficient rj at 
C — rj. The most relevant parameter (instead of rj) from 
the physical view point is the total number of quarks in 
the droplet Ng (as discussed above) and it is depicted 

to the left of each curve. (The variation of M, po and 
/tt could be considered as well instead of two mentioned 
parameters rj and C.) 
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TABLE I: Results of fitting by tlie Fermi distribution witli 
iV/ = 1, po (ch/fm^), Ro, t, ro (fm), b (fm-i), (MeV). 





Po 


-Ro 


b 


t 


ro 




V 


15 


0.34 


1.84 


0.51 


2.24 


0.74 


351 


0.65 


43 


0.43 


2.19 


0.52 


2.28 


0.75 


384 


0.73 


159 


0.46 


4.19 


0.52 


2.29 


0.77 


409 


0.78 


303 


0.47 


5.23 


0.52 


2.29 


0.78 


417 


0.795 


529 


0.47 


6.37 


0.52 


2.27 


0.79 


423 


0.805 


742 


0.47 


7.15 


0.52 


2.27 


0.79 


426 


0.81 



TABLE IL Results of fitting by the Fermi distribution with 
Nf = 2, po (ch/fm^), Ro, t, ro (fm), b (fm-^), m„ (MeV). 





Po 


Ro 


b 


t 


ro 


ma- 


V 


18 


0.81 


1.56 


0.37 


1.63 


0.57 


524 


0.7 


46 


0.9 


2.14 


0.37 


1.63 


0.6 


557 


0.75 


169 


0.93 


3.43 


0.36 


1.6 


0.62 


586 


0.79 


278 


0.94 


4.08 


0.36 


1.6 


0.62 


594 


0.8 


525 


0.94 


5.04 


0.36 


1.6 


0.62 


603 


0.81 


776 


0.94 


5.76 


0.36 


1.6 


0.63 


607 


0.815 



Analyzing the full spectrum of solutions obtained by 
scanning one can reveal a recurrent picture (at a certain 
scale) of kink-droplets which are easily parameterized by 
the total number of quarks Nq in a droplet and by the 
density po- These characteristics are evidently fixed at 
completing the calculations. The sign which allows us to 
single out these solutions is related to the value of droplet 
specific energy (see below). 

Table |T] exhibits the results of fitting the density p{r) 
with the Fermi distribution 



Ppir) 



Po 



1 



o{R.o-r)/b 



(34) 



where po is the density in origin, Rq is the mean size of 
droplet and the parameter b defines the thickness of sur- 
face layer t — 41n(3)6. Besides, the coefficient ro which 
is absorbed in the surface tension coefficient ([55]) . the a- 

1/3 

meson mass, Rq — r^Nq and the coefficient 77 at which 
all other values have been obtained are also presented in 
the Table H 

The curves plotted in the Fig. [10] and results of Table 
H] allows us to conclude that the density distributions at 
Nq > 50 are in full agreement with the corresponding 
data typical for the nuclear matter. The thicknesses of 
transition layers in both cases are also similar and the 
coefficient t-q with the factor 3^/^ included is in full cor- 
respondence with Fq. The values of cr-meson mass in 
TableUlook quite reasonable as well. However the corre- 
sponding quantities are strongly different at small quark 
numbers in the droplet. We know from the experiments 
that in the nuclear matter some increase of the nuclear 
density is observed. It becomes quite considerable for 
the Helium and is much larger than the standard nuclear 
density for the Hydrogen. 

Obviously, we understand the Thomas-Fermi approxi- 
mation which is used for estimating becomes hardly jus- 
tified at small number of quarks, and we should deal with 
the solutions of complete equation system (|28l) . However, 
one very encouraging hint comes from the chiral soliton 
model of nucleon [18], where it has been demonstrated 
tliat solving this system the good description of nu- 
cleon and A can be obtained. Then our original remark 
could be that the soliton solutions obtained in [Tsl] permit 
an interpretation as a 'confluence' of two kinks. Each of 
those kinks 'works' on the restoration of chiral symme- 
try since the scalar field approaches its zero value at the 



distance of ^ 0.5 fm from the kink center. Indeed, one 
branch of our solution corresponds the positive value of 
dynamical quark mass, and another branch presents the 
solution with negative dynamical quark mass (in three- 
dimensional picture the pseudo-scalar fields appear just 
as a phase of chiral rotation from positive to negative 
value of quark mass) . Such solutions develop the surface 
tension coefficient which is larger in factor two than the 
corresponding coefficient of single kink and as we believe 
signal some instability of a single kink solution. 

The similar results are obtained for two flavours Nf — 
2 (7 = 2NfNc = 12) assuming all dynamical quark 
masses of SU(2) flavour multiplet are equal. The solu- 
tions for the a-field and density distributions are similar 
to the corresponding results presented in Fig. [Tn]and Fig. 
[TTI The other data of fitting solutions are shown in the 
following Table [III As it is seen the characteristic ensem- 
ble density is approximately in factor two larger than 
the density of normal nuclear matter (remember again 
the factor 3). The characteristic values of cr-meson mass 
are slightly larger than for Nf = 1 and, consequently, 
the thickness of the transition layer is smaller almost in 
factor 1.4. The coefficient interrelating the mean size of 
droplet and the baryon (quark) number Fq ^ 0.8 is get- 
ting smaller. In principle, one can correct (increase) the 
surface layer thickness and the parameter Fq by decreas- 
ing the cr-meson mass but the ensemble density remains 
higher than the normal nuclear one. 



CD -3 




FIG. 12: The specific binding energy at Nf 
in MeV as a function of quark number Nq. 



1 and iVf 
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Fig. [n\ displays the specific binding energy of ensem- 
ble. It is defined by the expression similar to Eq. ([55)) 
in which the integration over the quark droplet volume 
is performed. The specific energy is normalized (com- 
pared) to the ensemble energy at the spatial infinity, i.e. 
in vacuum. Actually, Fig. [12] shows several curves in 
the upper part of plot which correspond the calculations 
with Nf — 1. The solid line is obtained by scanning over 
parameter rj and corresponds to the data presented in 
Table [Tj The dashed curve is calculated at fixed 77 = 0.4 

but by scanning over parameter M- It is clearly seen if 
the specific energy data are presented as a function of 
quark number Nq then the solutions, in which we are in- 
terested, rally in the local vicinity of the curve where the 
maximal binding energy - is reached. 

The similar solution scanning can be performed over 
the central density parameter po in origin. The corre- 

sponding data are dotted for a certain fixed M and pQ. 
It is interesting to notice that at scanning over any vari- 
able discussed a saturation property is observed and it 
looks like the minimum in Cb at Ng ^ 200-250. The 
results for the specific binding energy as a function of 
particle number are in the qualitative agreement with 
the corresponding experimental data. And one may say 
even about the quantitative agreement if the factor 3 (the 
energy necessary to remove one baryon) is taken into ac- 
count. Another interesting fact to be mentioned is that 
there exist the solutions of system ([28| with positive spe- 
cific energy. For example, for Nf = 2 such meta-stable 
solutions appear at sufficiently large 77 and with the den- 
sity parameter in origin equal po ~ P; = 0.157 ch/fm"^. In 
fact, the equation system p8[l represents an equation of 
balance for the current quarks circulating between liquid 
and gas phases. 

As a conclusion we would like to emphasize that in the 
present paper we have demonstrated how a phase tran- 
sition of liquid-gas kind (with the reasonable values of 
parameters) emerges in the NJL-type models. The con- 
structed quark ensemble displays some interesting fea- 
tures for the nuclear ground state (for example, an exis- 
tence of the state degenerate with the vacuum one) , and 
the results of our study are suggestive to speculate that 
the quark droplets could coexist in equilibrium with vac- 
uum under the normal conditions. These droplets man- 
ifest themselves as bearing a strong resemblance to the 
nuclear matter. Elaborating this idea in detail is a great 
challenge which will take a lot of special efforts and we 
do hope to undertake them in near future. 

Authors are deeply indebted K. A. Bugaev, R. N. Faus- 
tov, S. B. Gerasimov, E.-M. Ilgenfritz, K. G. Klimenko, 
E. A. Kuraev, A. V. Leonidov, V. A. Petrov, A. M. Sni- 
girev and many other colleagues for numerous fruitful 
discussions. 



Appendix A: Mean energy functional 

The free part of Hamiltonian 
Hq = — J dx q{x) {ijS/ + im) q{x) = 

= J-0j3\P^\ [cose {A+{p; s)A{p; s) - B{p; s)B+{p; s)) + 
+ sine {A+{-p; s)B+{p; s) + B{-p; s)A(p; s))] , 
contributes into the mean energy as 

Tr{eHo}- / 7^ \pi\ (l-cos0) + 

r dp (^1) 

+ J b4|cos6' [n{p)+n{p)] , 

where Hq ~ Hq/{V2Nc) is the specific energy. Natural 
regularization by subtracting the free Hamiltonian Hq 
contribution (without pairing quarks and anti-quarks) 
has been done in the first term of Eq. ([All) because 
in our particular situation this normalization in order to 
have the ensemble energy equal zero at the pairing angle 
equal zero turns out quite practical. It just explains a 
presence of unit in the term containing cos e. 

The Hamiltonian part responsible for in- 
teraction, qf^^f^qq'f^^^q' , provides four 
nontrivial contributions. The term 
Tt{pBB^ B' B'^} generates the following items: 
VUp, s)t^j€pVM T)tia^,Vsi{p, s) (the 
similar term but with the changes Q,T — >■ Q',T' which 
generates another primed quark current should be added) 
and -2F(Q,r) i°7^y(Q', T')F(Q', r')tVy(Q, T). 
Here (as in all other following expressions) we omitted 
all colour and spinor indices which are completely 
identical to those of previous matrix element. The 
term Tr{pBAA'~^ B'~^} generates the following nontrivial 
contributions: F(p, s)_<_"7^C/(q, t)I7(q, t)tS^^(P, «) 
-V{p, s)t"7^C/(P, S')C/(P, S')tVV'(p, s) 
-y(Q, T)t-rU{ci, t)U±q, t)t''rV{Q, T) 
+V(Q, r)r7'^C/(P, S')C/(P, S')tS''^(Q, T). Aver- 
aging Tr{pAA^ A' A'^} gives the contributions as: 
I7(P,5)i°7^C/(p,s)C7(p,s)tVC^(P,'S') (adding the 
similar term but with the changes P, S* — > P', S') 
and -2C7(P, SY'^-i^'UiV', 5')C7(P', S')t''YV{V, S). 
Another nontrivial contribution comes from av- 
eraging Tr{pA^ B^ B' A'} and it has the form 
V{q„T)t''YU{V,S)TJ{V,S)t''-^^'V{q„T). ah other 
diagonal matrix elements generated by the following 
terms Tr{pAA+P'P'+}, TT:{pBB'+ A'+ A'}, do not 
contribute at all (their contributions equal to zero). 
Similar to the calculation of matrix elements at zero 
temperature performed in Ref. Q we should carry 
out the integration over the Fermi sphere with the 
corresponding distribution functions in the quark and 

anti-quark momenta 10)^ ^ / (27r)-' ["(p) "-(-P)] 
deal with a finite temperature. All necessary formulae 
for the polarization matrices which contain the traces of 
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corresponding spinors could be found in Refs. and [s^. the result for mean energy density per one quark degree 
Bearing in mind this fact here we present immediately of freedom as 



\pi\ cos 0[n{p) + n{p)] + 2G f J^sm 



(27r)3 
dp 



(27r)3 
dq 



{6 - Om) Hp) + n{p)] 



J (27r)3 ' J (27r) 

I — 



(27r)3 

3 sin(0'-C) Hq)+n{q)] 1 + 



smit' — 



{e'-o'ji- 



(A2) 



(up to the constant unessential for our consideration 
here)[2^. It is quite practical to single out the colour fac- 
tor in the four-fermion coupling constant as G = 2G/Nc- 
Performing now the following transformations while in- 
tegrating in the interaction terms 

2 J dp f J dq- J dp f J dq r - J dp J dq^ 

= j dp f j dq{l- n- j dp{\-f) j dq, 

and changing the variables p o q in the last term we 



dq 

(2^ 



obtain 



j dp f j dq(l- n- j dp j dq (1 - /') 
- 1 dp (1 - /) I dq (1 - /') . 



Here the primed variables correspond to the momentum 
q. Then putting all the terms together we come to the 
equation fTT]) . 
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